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We derive a method to improve particle identification (PID) at high transverse 
momentum (pt) using the relativistic rise of the ionization energy loss (dE/dx) when 
charged particles traverse the Time Projection Chamber (TPC) at STAR. Electrons 
triggered and identified by the Barrel Electro-Magnetic Calorimeter (BEMC), pure 
protons(anti-protons) and pions from A (A), and K^ decays are used to obtain the 
dE/dx value and its width at given (3^ = p/m. We found that the deviation of the 
dE/dx from the Bichsel function can be up to 0.4cr (~ 3%) in p+p collisions at 
^/snn = 200 GeV taken and subsequently calibrated in year 2005. The deviation is 
approximately a function of (3^ independent of particle species and can be described 
with the function f{x) = A + „ ^ . The deviations obtained with this method are 
used in the data sample from p+p collision for physics analysis of identified hadron 
spectra and their correlations up to transverse momentum of 15 GeV/c. The ratio 
of e^/e^ (dominantly from 7-conversion) is also used to correct for the residual 
momentum distortion in the STAR TPC. 
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1 1 Introduction 

2 Identified hadron (tt^, K^, p{p)) spectra at higli pj- in p+p collisions provide a 

3 good test of perturbative quantum chromodynamics (pQCD) [1] and a baseline 

4 for studying color charge effect of parton energy loss in heavy ion collisions 

5 [2][3] . Hadron identification at high px was achieved with the ionization energy 

6 loss (dE/dx) at the relativistic rise using the STAR Time Projection Chamber 

7 (TPC). The dE/dx for a given charged-particle track is calculated using a well- 
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8 known method called "truncated mean". The dE/dx values of hits associated 

9 with the track have a typical Landau tail. The hits with the top 30% of high 

10 dE/dx values are discarded and an average of dE/dx value from the rest 

11 of the hits is derived for that track [9]. The dE/dx for a given particle at 

12 low momentum decreases with increasing momentum to reach a minimum 

13 ionization, then increases due to the relativistic rise. For a minimum ionizing 

14 particle (MIP) the dE/dx resolution in the STAR TPC is 6-8% for a track 

15 with the maximum of 45 sampled dE/dx points. The pions are well separated 

16 from the rest of the particles {e,K,p) at p^ of 0.3 to 0.6 GeV/c. This clear 

17 separation has been used to calibrate the TPC dE/dx without other means of 

18 identification. It provides the fixed points for extrapolating the dE/dx function 

19 to higher momentum. In the thin material (TPC gas) the Bichsel function has 

20 proved to be a very good approximation for the dE/dx curves and has been 

21 adopted by STAR as a standard method of predicting the dE/dx value for 

22 charged hadrons in all momentum ranges [1]. 

23 The Landau function is an approximation which does not include features 

24 related to atomic structure. With the Bichsel function |1], a decrease of the 

25 relativistic rise of dE/dx with increasing segment length x is seen and param- 

26 eterized empirically. This effect, gas multiplication gains and noise of TPC 

27 electronics and pileup in high luminosity environment may make the dE/dx 

28 deviate from the Bichsel function. In the relativistic rise region, the dE/dx sep- 

29 arations among vr^, K^ and p{p) are about l-3cr with the dE/dx amplitude of 

30 pions the highest and that of protons the lowest. Pions are the dominant par- 

31 tide species for inclusive and jet hadrons, and they shadow kaons and protons 

32 in the dE/dx distribution. In a given momentum slice, clear peak separations 

33 of these three hadrons are not possible. This results in large systematic er- 



34 rors due to the uncertainty of the dE/dx positions. From the minimum-bias 

35 triggered p+p colhsions in the year 2003, the p(p) spectra were measured to 

36 Pr < 7 GeV/c with significant systematic errors due to the uncertainties in 

37 determining the mean dE/dx position for proton and kaon [3]. Knowledge of 

38 the precise dE/dx positions for those hadrons is important to understand the 

39 efficiencies of PID selection and to reduce the systematic uncertainty in iden- 

40 tified hadron yields. In order to improve the particle identification at high pt, 

41 we develop a method to locate the dE/dx positions for different hadrons with 

42 good precision. 

43 The p+p collision events with enhanced high-p^ charged particles used in this 

44 analysis were obtained from online jet trigger by the BEMC with < rj < 

45 1 in the year 2005, and full azimuthal coverage. A -^ p + n^ (A— ^p + 7r+) 

46 and Kg -^ tt+tt^ are reconstructed by their decay topology to identify their 

47 decay daughters - charged pions and protons. The identified electrons, pions 

48 and protons provide the necessary distinct dE/dx positions and widths as 

49 function of /37. The deviation from the prediction of the Bichsel function was 

50 used to correct for the dE/dx fit when extracting pion and proton yields from 

51 the charged hadrons in an inclusive hadron distribution or in a jet. The same 

52 method can be applied to p+p, d+Au and A+A collisions in STAR. 

53 Physics goals of the STAR Experiment at RHIC in recent (and future) years [5] 

54 drive the need to operate the STAR TPC at ever higher luminosity, and in- 

55 crease the ionization density in the TPC gas [6J. The resulting ionic space 

56 charge and grid leakage introduce field distortions in the detector, which sys- 

57 tematically shift the reconstructed momentum of positive and negative parti- 

58 cles in oppositive directions. The effect is expected to grow as a function of 

59 Pt- star has developed a method for correcting the track reconstruction due 



60 to space charge distortion. Performance of the corrections can be assessed by 

61 examining the distribution of signed DCA (Distance of Closest Approach of 

62 a primary track to the colhsion vertex) as a function of luminosity. In this 

63 paper, we use electrons and positrons dominantly from gamma conversion to 

64 check and correct for potential residuals from the distortion. We expect e~ je^ 

65 to be unity independent of p^ since 7 — >■ e~ + e+ and a significant fraction 

66 of leptons from heavy- flavor decays are also expected to be close to unity [7] . 

67 The high-statistics data set from BEMC trigger is ideal for such a study. 



68 2 Experimental Setup and Data Analysis 

69 2. 1 Detectors and Datasets 

70 The Solenoidal Tracker at RHIC (STAR) [8] is a powerful detector with a large 

71 and uniform acceptance capable of tracking charged particle and providing 

72 particle identification (PID). Its main tracking detector - TPC [9] covers full 

73 azimuthal angle ( A0 = 2ti ) and -1.3 < r/ < 1.3 in pseudo-rapidity, and 

74 provides identification of the charged hadrons by measuring momentum and 

75 dE/dx separation of the charged particles. The STAR TPC has 45 rows of 

76 pads with the first pad row at a radius of 60 cm and last at 180 cm, providing 

77 a maximum of 45 sampled dE/dx points with 78 cm track length in a PIO 

78 gas (90% Argon and 10% CH4) while Track lengths of up to 110 cm are 

79 possible. In order to obtain the dE/dx position of protons at the relativistic 

80 rise range, protons (anti-protons) are selected from A (A) decays. Similarly, 

81 decay pions from Kg can be used to get dE/dx and dE/dx positions of pions 

82 for 0.2 < pt < 3 GeV/c. The precise track geometry allows the particle 



83 trajectories to be extrapolated to form a VO away from the collision vertex [10]. 

84 This provides dE/dx information from samples of pure protons and charged 

85 pions. 

86 In addition, the sub-detector BEMC [IT] is also used to enhance electron and 

87 positron yields and to confirm their identification. This study used about 16.7 

88 million p+p collision events recorded through online High Tower (HT) and 

89 Jet Patch (JP) triggers in the BEMC. The HT trigger condition required the 

90 energy of a single calorimeter tower (A?7 x A0 = 0.05 x 0.05) to be at least 

91 2.6 (HTl) or 3.6 (HT2) GeV, and JP trigger required the total energy of one 

92 patch of towers (Ar/ x A0 = 1.0 x 1.0) to exceed 4.5 (JPl) or 6.5 (JP2) 

93 GeV [12] . 

94 2.2 dE/dx distribution in the TPC at high px 

95 To formulate the dE/dx ^ distribution and its associated Bichsel function for 

96 PID, we need to define a few terms. The normalized dE/dx is defined as 

97 nax = ^og[{dE/dx)Y/Bx]/o-x, where X,Y can represent vr^, K^, p(j>) or e^; 

98 Bx is the expected mean dE/dx of particle X; and ax is the \ia{dE/dx/Bx) 

99 resolution of the TPC p^l4fl5] . All the quantities are calculated in a track- 

100 by-track basis. Fig. [1] shows na.„ distribution of all the charged particles for 

101 3.75<]9t<4.0 GeV/c and 8.0<pt<10.0 GeV/c at | r^ |< 0.5. A sum of eight 

102 Gaussian functions is used to fit this distribution with thirteen parameters 

103 in order to obtain the identified hadron yields. Each Gaussian describes one 

104 dE/dx distribution for a charged particle species. The thirteen parameters are 

105 peak positions relative to pion peak ( naj, ncr^-naj, ncr^-ncr^, na^-na^ ) , eight 



^ dE/dx is used to represent the "track descriptor" C defined on p. 170 of [1] 



106 yields for the charged particles and one common Gaussian width. With ideal 

107 calibration, na'^ should be a normal Gaussian distribution centered at zero and 

108 with a width of unity. The normality of Gaussians depends on the precision of 

109 ax in a track-by-track basis. As presented in Ref. [3] (section 6.2 and Fig. 21, 
no 22, and 23), the ax depends on the track length (t) in the dE/dx sample. An 

111 empirical parameterization yields a^ = ctqA* where 0.45 < s < 0.55 [1]. Fig. [2] 

112 shows the dE/dx resolution as a function of track length for this specific run 

113 in p+p collisions at center of mass energy of 200 GeV in year 2005. The red 

114 solid curve is a power-law fit with s = 0.52 in the range of 40 < t < 80 cm, 

115 where the majority of the tracks come from. A polynomial function was also 

116 performed at larger range and was shown as dashed line. 



117 Fig. [T] shows that the pion dE/dx position is different from the Bichsel func- 

118 tion. This means that the dE/dx calibration is not perfect and it also implies 

119 that dE/dx position of other particles relative to that of the pions may be off 

120 from their expected values. In order to improve the particle identification and 

121 reduce the systematic uncertainty in identified particle yields obtained from 

122 dE/dx [2||3] . we study in details the precise dE/dx positions of all charged par- 

123 tides using the enhanced electron by the BEMC, pure proton decayed from 

124 from A and pion decayed from Kg in the TPC. Once all the dE/dx positions 

125 and widths for all the charged hadrons are obtained by other means, we will be 

126 able to constrain better the Gaussian fits, and understand the efficiency and 

127 contamination better in the case of PID selections for other physics analysis. 



128 2.3 Electron identification by the BEMC 

129 Although electron dE/dx is relatively far away from those of the other charged 

130 particles, the electron yields are orders of magnitudes lower than the yields 

131 for pions. In order to identify the electron and obtain its dE/dx position, a 

132 dataset with a special trigger based on the energy deposited in the BEMC 

133 tower is used to enhance the yield of electrons relative to other particles. Ad- 

134 ditional hadron rejection is achieved from the shower shape and position from 

135 a Shower Maximum Detector (SMD) [TT], which is located at ~5.6 radiation 

136 lengths depth in the BEMC. A unique feature of the SMD is its double layer 

137 design which makes it possible to reconstruct the shower as two-dimensional 

138 image, so that it can provide fine spatial resolution in and rj directions and 

139 reject hadrons according to the different shower shape between hadrons and 

140 electrons. We require p/E to be 0.3< p/E <1.5 where p is track momentum in 

141 the TPC and E is the energy of the BEMC tower, and the shower hits of the 

142 SMD in rj and (p direction to be ra^ > 2 and (p direction n^ > 2, respectively. 

143 The different positions between hit and track projection in (p and z directions 

144 are restricted to be | (pdist |< 0.01 rad and | Zdist |< 2 cm [TT]. The na^ dis- 

145 tribution for track after these cuts are shown on Fig. [3] for 3.75 < Pt< 4.0 

146 GeV/c and at much higher p^ range. With about 1.5-3o" separation between 

147 electron (positron) and other hadrons, electron position and yields could be 

148 obtained from the eight-Gaussian function as above. From the same fit, we 

149 can also obtain the dE/dx positions for pions, which are labeled as "pion with 

150 EMC" for later use. 

151 To correct for and to assess the systematic errors due to the residual mo- 

152 mentum distortion, we use e~/e~^ ratio as a function oi Pt- Fig. H] shows the 



153 ratio as a function of px obtained from the BEMC triggered data described as 

154 above. If the electron and positron yields are a modified power-law function 

155 ifipx) oc {po + Pt)^^) without any distortion, the distortion due to space 

156 charge in the TPC shifts all negatively charged tracks from pt value to a 

157 higher px + A x p\ while it shifts all positively charged tracks from p^ value 

158 to a lower px — A x p]. \Q\. The correction form is according to the effect of 

159 sagitta displacement due to space charge on momentum distortion [16]. Data 

160 points are fitted by the following function f{pT) = ( 2'67+^^-Al^^ )^^'"^' where 

161 A* p\ means Ap^ affected by charge distortion, and po = 2.67 and n = 11.4 

162 are parameters obtained from the inclusive electron spectra [7j. The px depen- 
ds dence of ratios indicates that the momenta of the charged particles obtained 

164 from the TPC tracking still systematically shifted away from their true value 

165 due to the space charge distortion. We note that the obtained distortion char- 

166 acterized by parameter A is only about 2a from zero and this results in about 

167 1.3% momentum shift for a single particle track at pt = 15 GeV/c. 



168 2.4 Proton and pion from VO reconstructed in the TPC 

169 In order to obtain the dE/dx position of proton {na^) for (3"^ > 4, protons 

170 (anti-protons) are selected from A (A) through A -^ p + n^ (A^p + 7r+) 

171 decays, because it is difficult to get dE/dx position of proton by h'^ — h~ [13] 

172 with this data sample due to low statistics and small difference of yields of 

173 proton and anti-proton. At the same time, pions decayed from K^ through 

174 Kg — i> vr"^ + 7r^ decay can be used to get dE/dx and dE/dx positions of pion 

175 (nal) at 0.2<pt<S GeV/c. First, Kg and A are selected by topological cuts 

176 on a secondary vertex ^U\ according to long decay length of Kg (cr=2.6 cm) 



177 and A (cr=7.89 cm). Fig. [5] shows the invariant mass distribution for K^ (top 

178 panel) and A (bottom panel). Pure K^ and A (A) are selected via their in- 

179 variant mass cuts, 0.485<M(K° )<0.505 GeV/c^ and 1.112< M(A) < 1.12 
iBo GeV/c^, and their daughter particles (vr,p(p)) with high purity are obtained 

181 to derive n&^ and na^. Fig. [6] shows na^ distribution of pions decayed from 

182 K^ (upper panel) fitted by a Gaussian function and protons decayed from A 

183 (lower panel) fitted by a 2-Gaussian function. The protons from A decay have 

184 higher background (signal-to-background ratio = 9:1) and a second Gaussian 

185 representing the pion contamination is necessary. Meanwhile, the pt depen- 

186 dence of the dE/dx width for protons and pions from the fits are shown in 

187 Fig. [7l which also shows dE/dx width for electron. The width is consistently 

188 smaller than unity (0.868 ± 0.004). This means that the dE/dx resolution is 

189 about 13% better than the prediction and the separations among particles 

190 are better than what we expected. This may be due to the run to run vari- 

191 ation (luminosity, beam background, etc.), and our calibration only sampled 

192 a small fraction of data(~ 5%). The open triangles in Fig. [7] shows the sole 

193 effect of the smearing of the dE/dx peak positions between proton and pion 

194 due to the variation of the track quality (cr^). This effect (< 0.2o"), compared 

195 quadratically to the dE/dx width, contributes < 4% to the final dE/dx width. 



196 2.5 dE/dx deviation vs Pj 

197 With identified pion, proton and electron mentioned above, the experimental 

198 results on the deviation of the normalized dE/dx {na^) relative to the Bichsel 

199 theoretical values, as a function of /37 are shown on Fig. [HI Since there is 

200 almost no particle species dependence of dE/dx, we can describe it with a 

10 



201 function of f{x) = A + qt:^- The fit parameters are listed in Tab. [T] With 

202 this, we can determine the dE/dx positions and widths of any given charged 

203 particles to better than < O.lo" or < 1%. 

204 There are two ways to correct for this effect in the data. One can attempt 

205 to understand the origin of this deviation and correct for the effect at the 

206 hit level (the amplitudes of the ionization signal in each pad and row). This 

207 requires re-processing of the hits and reconstructing tracks from scratch. In 

208 order to take advantage of the existing compressed dataset with tracking in- 

209 formation only, we apply the corrections to the dE/dx Gaussian function for 

210 each particle species without modifying the dE/dx itself. This empirical after- 

211 burner is applied in each px and rapidity bin to directly extract particle yields 

212 required by the physics analyses. Fig. M shows px/'ma.ss dependence of the 

213 normalized dE/dx deviation, which is fitted by the same function as for the 

214 case of p/mass. The parameters from the function are shown in Tab. [TJ With 

215 the corrected deviation, differences of dE/dx between pion and other charged 

216 particles {na^ — na^, na^ — na'^ and na'^ — na'^) are calculated and compared 

217 with theoretical values and shown in Fig. [TOl Clear offsets are depicted in 

218 Fig.[IOl 

219 3 Summary 

220 The dE/dx positions and widths of charged particles have been precisely de- 

221 termined using BEMC triggered events with enhanced electron content and 

222 pure samples of pions and protons from KOs and Lambda decays respectively. 

223 Their deviations relative to theoretical values versus /37 and pT/'ma.ss are de- 

224 scribed by the function f{x) = A+ qt^ very well. With this method, dE/dx 
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225 positions of charged particles are re-calibrated to be better than O.lcr. The 

226 particle identification of charged hadrons is thus improved, and the uncer- 

227 tainty is reduced significantly. These are important steps toward fulfilling the 

228 physics goals of the STAR experiment at RHIC in the future. 
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Fig. 1. ncr^ distributions at 3.75<pt<4.0 GeV/c and 8.0<pt<10.0 GeV/c for the 
positive (left panel) and negative (right panel) particles. The solid lines are the 
8-Gaussian function, which is a sum of the invidivual Gaussian functions from pion 
(dashed line), kaon (dot-dashed line), proton (dotted line), and electron (long-dash 
line). The solid vertical lines are the extracted pion dE/dx positions, while the 
dashed vertical lines are the previously expected positions. 



Resolution versus Track Length 




Track Length (cm) 



Fig. 2. dE/dx resolution {ct-k) as a function of track length. The red-solid line is a 
uq/I^-^'^ power-law fit, and the dashed line is a polynomail function upto 4th power 
in In (t). 
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Fig. 3. na^ distribution with high-tower trigger at 3.75<pT<4.0 GeV/c (left panel) 
and 8.0<pt'<10.0 GeV/c (right panel). Also shown is the result from an 8-Gaussian 
fit. Both electron and positron yields are enhanced significantly by the trigger. 



2.5 



1.5 






1 - 



0.5 



- 






-_ 








" r-*+- 


«±J 


ry-H 








T 






" , , 1 





4 6 8 10 12 14 

p^(GeV/c) 



Fig. 4. The ratio of e /e^ as a function oipr- The curve is a power-law fit described 
in the text. 

Table 1 

Fit parameters from momentum andpr dependence of ncr^ with the function f{x) = 



A + 



B 



C+x^ 



parameters 


xVndf 


A 


B 


C 


p dependence 


95/49 


-0.423 ± 0.015 


235 ± 23 


464 ± 37 


Pt dependence 


94/48 


-0.443 ± 0.015 


234 ± 23 


444 ± 35 
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Fig. 5. Invariant mass distributions of i^^ Fig. 6. na!^ distributions of pions from 

(upper panel) and A (lower panel) with K^ (upper) and protons (lower) from A. 

at least one daughter particle (vr or pro- The solid line is for pions, and the dashed 

ton) at 3.0<pt<3.25 GeV/c. one for protons. 
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Fig. 7. "Pt dependence of dE/dx na^ Fig. 8. The dE/dx deviation as a func- 

width for proton from A decay (crosses), tion of (5^. The crosses are for protons 

for pion from K^ decay (circles), and from A, the filled dots are for pions from 

for electron enhanced by the EMC (di- K^ , and open circles and stars are for 

amond). The open triangle is the smear- pions and electrons from HT trigger, re- 

ing between proton and pion dE/dx peak spectively. 
positions due to the variation of the track 
quality (o-^). 
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Fig. 9. The dE/dx deviation as function Fig. 10. The relative dE/dx peak posi- 
of pT/mass (similar to Fig. [5]). tions of na!^ , na^, na'^ as a function of 

Pt- The solid dots depict theoretical val- 
ues, and open ones are re-calibrated re- 
sults. 



